function[x1,y1,z1] = ga-gm(x,y,z),

theta = (90-79.30)*pi/180; %(lat)
phi = (288.59)*pi/180; % (log)

comatrix = [cos(theta)*cos(phi), cos(theta)*sin(phi), -sin(theta);...
           -sin(phi),           cos(phi)                       0;...
           sin(theta)*cos(phi) sin(theta)*sin(phi)  cos(theta)];
   
   
   